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array  periodicity  (Bloch  oscillator) . 
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2. 


The  steady  increase  of  large  scale  integration  requires 
smaller  and  smaller  structures  with  thinner  and  thinner 
layers.  Simultaneously  the  number  of  defects  has  to  be 
decreased  so  that  any  kind  of  defect  analysis  is  very 
important  for  improvement  of  the  production  of  devices.  One 
key  issue  is  the  quality  of  the  Si/Si02  interface.  Our  own 
results  of  the  recent  years  have  demonstrated,  that  defects 
at  the  interface  are  identified  and  quantitatively  measured 
by  spot  profile  analysis  of  LEED.  Whereas  TEM  and  XPS 
provide  also  some  qualitative  and  only  limited  quantitative 
results,  only  LEED  was  so  far  able  so  measure  systematic 
variations  of  atomic  roughness  at  the  interface  due  to 
variation  of  oxidation  parameters  (see  our  final  report 
DA  JA  37-81-C-0054  and  references  /!/  and  /2/) .  Additionally 
a  direct  correlation  of  atomic  roughness  to  electronic 
properties  like  mobility  and  interface  state  density  could 
be  shown. 

The  purpose  of  the  present  investigation  has  been,  to 
analyze  the  atomic  roughness  at  the  interface  in  more 
detail,  so  that  the  quantitative  description  is  more 
complete.  The  investigations  should  be  expanded  from 
Si  (111)  to  the  technologically  more  important  Si (100)  which 
is  more  complicated  due  to  the  equidistant  sublattices  of 
the  diamond  lattice  in  the  [100]  direction.  With  the  better 
structural  information  a  more  detailed  correlation  to 
electrical  data  like  mobility,  ESR  and  radiation  damage 
should  be  possible.  The  maximum  field  strength  for 
electrical  breakdown  should  be  related  to  atomic  roughness 
too . 

The  results  have  been  obtained  in  part  by  LEED 
investigations  in  our  laboratory.  Those  studies  have  been 
very  successful  and  have  a  lot  of  results  not  even  expected 
in  the  proposal.  For  the  comparative  studies  we  needed 
outside  nclp  for  obtaining  oxidation  and  electrical 
measurements  with  commercial  standard.  Those  measurements 
were  badly  delayed  for  reasons  we  could  not  influence,  so 
that  this  part  of  the  proposal  is  only  fulfilled  in  part,  it 
is  to  be  continued,  until  the  missing  measurements  are 
completely  added. 

The  results  will  be  described  according  to  the  main  aspect 
of  each  special  study. 


Development  of  Roughness  Analysis 


Since  in  this  study  for  the  first  time  the  newly  developed 
high  resolution  LEED  system  has  been  used,  it  was  necessary 
to  develop  the  evaluation  procedures,  so  that  the  improved 
set  of  data  is  used  not  only  for  a  higher  accuracy  (a 
factor  of  up  to  20  has  been  obtained) .  It  has  been  rather 
desirable  to  find  also  new  kind  of  informations  not  yet 
available  so  far. 
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3.1  Terrace  width  probability  distribution 

Up  to  now  only  fitting  procedures  have  been  used  to 
determine  the  probability  distribution.  Here  just  a 
distribution  had  to  be  assumed,  then  the  autocorrelation  and 
its  Fourier  transform  was  constructed  and  compared  with  a 
measured  profile  (see  ref.  111).  So  by  trial  and  error  a 
best  fit  was  approximated.  A  straightforward  derivation  has 
been  found  in  ref.  /4/,  which  provides  a  direct  transform  of 
the  spot  profile  into  a  probability  distribution.  The 
details  are  found  in  ref.  /4/.  As  a  first  success  for  the 
first  time  it  could  be  shown,  that  the  contribution  of  the 
background  to  the  profile  is  of  decisive  influence  onto  the 
quantitative  evaluation  (ref.  / 4 / ) . 


3.2  Vertical  structure 

So  far  only  the  lateral  structure  could  be  derived  out  of 
the  spot  profile  .  The  high  resolution  profiles,  however, 
showed,  that  in  most  cases  the  profile  consists  out  of  a 
central  sharp  spike  and  a  well  separated  broad  shoulder  /6/. 
It  could  be  demonstrated,  that  the  fraction  of  intensity  in 
the  spike  out  of  the  integrated  spot  intensity  varies 
characteristically  with  energy.  This  variation  is  easily 
described  within  the  kinematic  approximation  and  contains 
just  the  vertical  probabilities  (the  probability  of  finding 
a  surface  atom  at  a  given  vertical  level)  .  Consequently  it 
is  now  for  the  first  time  possible  to  derive  those 
probabilities  (at  least  their  autocorrelation)  161.  If  a 
small  number  of  levels  is  involved  or  if  a  symmetric 
distribution  is  assumed,  a  unique  evaluation  is  possible, 
which  enables  an  excellent  atomic  view  of  surface  structure. 
The  details  of  analysis,  its  experimental  proof  and  its 
application  to  MBE  is  found  in  ref.  /6/. 


3.3  RMS-Deviation  from  average  surface 

If  all  vertical  probabilities  pi  are  known,  the  rms- 
roughness  _  is  easily  derived 

l1  =  I  piiJ  -  Z  ( pi  •  i ) 2 


Due  to  ambiguity  for  pi  data  at  a  large  number  of  layers  it 
is  desirable  to  have  a  direct  and  model  independant 
determination  of  the  rms-roughness .  Such  a  determination  is 
possible  with  the  energy  dependence  of  the  function  G(R  ), 
which  is  the  above  function  to  describe  the  fraction  of  the 
integral  peak  intensity  within  the  central  spike.  If  the 
function  GOO  is  represented  by  its  Taylor  expansion  at  an 
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in-phase-condition,  already  the  second  coefficient  provides 
the  roughness  i. 


G  (K  )  a  1  -  ^2K  2 


Since  the  experiment  yields  directly  the  roughness  out  of 
the  curvature,  the  determination  is  completely  independent 
of  any  model  or  assumption  /9/.  Additionally  the  roughness 
±  and  the  correlation  L  are  the  quantities  needed  to 
describe  the  mobility  due  to  roughness  scattering 
quantitatively.  Since  both  quantities  are  directly  derived 
out  of  the  LEED  measurement,  the  comparison  of  mobility  and 
roughness  is  now  possible  without  fitting  parameter. 


3.4  Chemical  inhomogenities 

It  is  well  known,  that  any  defect,  structural  or  chemical, 
may  produce  a  broadening  of  a  spot  profile.  The  defects  due 
to  steps  are  easily  separated  by  the  periodic  variation  of 
the  profile  with  energy  { / 1  —  3 / )  .  It  is  desirable,  to  know 
simultaneously  any  other  defects  at  the  surface.  A  general 
concept  has  been  recently  developed  with  the  inclusion  of 
two  different  formfactors,  which  besides  steps  also  gives 
rise  to  a  shoulder  (completely  described  in  the  diploma 
thesis  of  J.  Falta  (1988) ,  to  be  published  soon) .  This 
procedure  has  been  used  to  distinguish  at  the  Si02/Si 
interface  between  steps  and  other  defects,  which  will  be 
discussed  in  the  next  chapter  (see  also  /9/)  . 


4.  The  Annealing  of  the  Si0a/5i(lll)  interface 


4.1  Experimental  results 

Commercially  polished  Si  (111)  wafers  have  been  oxidized  in 
dry  oxygen  at  800° C  until  a  10  nm  thick  oxide  has  been 
grown.  Some  wafers  have  been  annealed  in  nitrogen  at  800° C 
or  1000° C  for  different  periods  of  times.  The  purpose  has 
been  to  determine  defects  with  high  resolution  and  to 
measure  changes  due  to  annealing. 

As  developed  in  earlier  studies  /l/  the  oxide  has  been 
removed  in  HF  and  transferred  under  a  droplet  of  methanol 
into  the  ultrahigh  vacuum,  so  that  a  reoxidation  in  air 
could  be  avoided.  The  high  resolution  SPA-LEED  was  similar 
to  that  described  in  /6/.  Due  to  the  high  transfer  width  of 
more  than  100  nm  (at  least  5  times  better  than  in  earlier 
studies  of  ref.  /I/  and  to  the  high  dynamic  range  of  5 
orders  of  magnitude  quite  a  few  informations  have  been 
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obtained,  which  are  qualitatively  and  quantitatively  far 
beyond  the  earlier  studies. 

The  spot  profiles  have  been  measured  along  the  reciprocal 
lattice (00) -rod  between  the  555  and  666  Bragg  spots,  that  is 
between  90  and  140  eV,  so  that  the  phase  varies  between 
5  <  S  <  6  with  2n.S  the  phase  shift  of  scattering  at  ad¬ 
jacent  levels.  The  spot  profile  at  the  out-of-phase 
condition  (S  =  5.5!  is  shown  in  fig.  4.1  for  the  not- 
annealed  wafer.  For  comparison  purposes  also  as  a  solid  line 
the  profile  is  shown  for  a  well  annealed  surface,  which 
reproduces  the  instrumental  resolution.  The  difference  of 
the  two  curves,  which  after  two-dimensional  integration  is 
much  higher  than  the  central  spike,  shows  the  broadening  due 
to  defects.  For  quantitative  evaluation  the  measured  profile 
is  fitted  by  a  sum  out  of  four  components:  the  central  spike 
(given  by  a  6-function  convoluted  with  the  instrument 
response  function) ,  two  Lorentzian  profiles  with  different 
half  widths  and  a  constant  background. 

The  fou’'  contributions  are  found  by  a  logarithmic  plot  of 
the  out-of-phase  pair  correlation,  which  is  the  fourier 
transform  of  the  measured  profile.  Whereas  background  and 
central  spike  are  separated  by  the  pair  correlation  value  at 
zero  and  the  values  at  high  distances  (more  than  some 
hundred  atomic  distances) ,  the  two  Lorentzians  are  now 
available  as  linear  portions  of  the  pair  correlation. 
Whereas  the  constant  background  is  not  further  discussed  as 
being  due  to  thermal  diffuse  scattering  or  random 
distributed  point  defects,  the  other  three  ones  are  studied 
carefully  by  measurement  at  many  different  energies.  Fig. 
4.2  shows  the  spot  profile  of  the  same  sample  as  in  fig. 
4.1  here  for  in-phase-condition  (S  =  5.0).  It  is  easily 
seen,  that  the  profile  consists  not  just  of  a  central 
spike  (and  a  constant  background) ,  one  broad  Lorentzian  has 
to  be  added.  By  comparison  with  fig.  4.1  it  is  found,  that 
the  same  broad  shoulder  persists,  only  the  small  Lorentzian 
shoulder  of  fig.  4.1  is  gone.  It  is  therefore  concluded  that 
the  small  shoulder  is  due  to  steps  and  the  broad  shoulder 
due  to  other  defects.  For  a  clear  distinction  and  quantita¬ 
tive  evaluation  the  integral  contributions  of  the  three 
components  (neglecting  the  background)  are  measured  relative 
to  their  total  intensity.  The  results  are  shown  in  fig.  4.3 
with  the  contribution  Is  of  the  narrow  Lorentzian  shoulder 
and  Ib  of  the  broad  Lorentzian  shoulder.  The  shape  has  been 
except  scaling  identical  for  all  samples,  so  that  an 
evaluation  is  consistently  possible.  The  broad  shoulder  does 
not  disappear  for  in-phase-condition,  it  is  therefore  due  to 
defects,  described  by  a  different  form  factor  fz  (as 
explained  in  chapter  3) .  It  is  described  for  most  energies 
within  the  approximation  of  a  very  small  form  factor  of  the 
defect  areas,  so  that  a  constant  value  for  Ib  is  obtained 
for  most  of  the  energies.  Only  close  to  S  =  5.5  larger 
values  are  obtained.  Here  the  intensity  for  the  clean 
silicon  is  very  low,  so  that  the  approximation  fi  >>  f2  does 
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no  more  hold.  The  deviation  is  therefore  understood.  The 
constant  value  is  used  to  derive  the  coverage  of  the  surface 
with  defect  areas  of  low  form  factor,  because  its  value  is 
directy  the  coverage  6b  with  different  formfactor  fs . 

The  narrow  shoulder  shows  a  nearly  sinusoidal  variation  with 
energy,  which  is  characteristic  for  a  two-level  stepped 
surface.  For  a  quantitative  evaluation  not  its  contribution 
to  the  total  integrated  intensity  is  useful,  according  to 
the  theory,  rather  its  contribution  to  the  sum  out  of 
central  spike  and  narrow  shoulder  should  be  used.  Fig.  4.  4 
shows  its  complement,  that  is  the  integral  of  the  central 
spike  over  the  sum  out  of  this  spike  and  the  narrow 
shoulder,  which  reproduces  the  sinusoidal  shape  to  a  more 
large  extent.  This  function  is  evaluated  according  ref.  8 
as  a  perfect  two-level-system  with  a  coverage  in  one  level 
of  about  35  %.  {see  table  1) .  The  rms  roughness  as  explained 
in  chapter  3  is  also  shown  in  table  1. 

The  lateral  distribution  of  both  steps  and  the  other  defects 
is  derived  out  of  the  shape  and  width  of  the  profiles.  Since 
both  profiles  have  exactly  a  Lorentzian  shape,  their  pair 
correlation  is  given  by  an  exponential  function.  Therefore 
the  length  or  size  distribution  is  given  by  the  geometric 
distribution 


1 

P  (D  =  -  (  1 

<r> 


l 

<r> 


r  -  l 

) 


where  P  (D  is  the  probability  of  finding  a  terrace  or 
domain  of  width  r  (in  multiples  of  the  lattice  constant  a) . 
Just  one  parameter,  the  average  width  <r>  is  sufficient  for 
description  of  the  length  distribution  for  each  entity 
(lower  or  upper  level,  defect  free  or  distorted  areas).  The 
shoulder  in  a  one  dimensional  description  is  given  by  the 
function 


1-5  4  9  (1  -  ©)  sin2  (K/d/2) 

Go  {  K ...  K  |  )  —  1  ’  *  ^ - 

1  +  5  1  +  S!  -  2  5  cos  ( K  ||  a) 


which  may  be  described  by  a  sum  of  identical  Lorentzians  at 
each  position  of  a  diffraction  spot.  The  parameter  5 
describes  the  expontential  decrease  of  the  pair  correlation 
*(n)  -  exp(-n/j>)  via  a  correlation  length 
a 


y 


In  5 


The  coverage  ©  as  derived  above,  and  the  parameter  8  provide 
the  average  width  <  To  >  and  <  Ti  >  for  the  two  levels  or 
the  two  kind  of  surface  patches 


<  To  >  =  (©  (1  -  (5)  )-* 

<  n  >  =  (  (1  -  ©)  (1  -  3 ) )  - 1 


Those  numbers  are  shown  in  table  4.1  and  4.2. 


4.2  Discussion 

The  results  summarized  in  table  4.1  and  4.2  show,  that  the 
interface  is  extremely  smooth.  The  roughness  is  given  by 
about  half  the  layer  distance.  Only  a  change  between  two 
levels  along  the  interface  is  observed.  The  average  terrace 
width,  as  derived  out  of  the  energy  dependent  shoulder  is 
fairly  large.  Both  the  roughness  and  the  terrace  width  are 
nearly  not  affected  by  heat  treatment. 


An  independent  and  completely  new  information  is  derived  out 
of  the  broad  shoulder  with  the  nearly  energy  independent 
amplitude.  It  is  concluded  that  at  the  surface  only  part  of 
the  surface  consists  of  well  ordered,  clean  silicon  (the 
percentage  is  given  in  table  2,  as  derived  out  of  the 
intensity  of  the  broad  shoulder  at  energies  close  to  the  in- 
phase  condition) .  The  rest  of  the  surface  is  covered  with 
domains  of  different,  for  most  energies  much  lower  form 
factor.  Those  domains  again  have  a  size  and  distance 
distribution  as  given  by  the  geometric  distribution.  There¬ 
fore  their  average  length  can  be  given  in  table  2. 

Possible  origins  for  those  domains  of  low  scattering  power 
(or  form  factor)  are  discussed  now.  It  may  be  due  to  areas 
with  different  chemical  composition.  It  has  been  shown,  that 
a  fraction  of  a  monolayer  of  oxygen  and  carbon  is  found  at 
those  surface.  On  the  other  hand  structural  disorder  in  some 
areas  may  work  the  same  way.  It  may  Le  a  combination  of  the 
two,  that  oxygen,  embedded  a  little  bit  deeper  during 
oxidation  produces  some  structural  disorder.  The  large  half 
width  of  the  shoulder  means,  that  those  distorted  areas  are 
fairly  small,  only  3  to  8  atomic  distances  in  average.  The 
distribution  again  is  geometric,  pointing  to  a  random 
process  of  generation  without  correlation.  In  big  contrast 
to  the  step  structure  the  average  domain  size  and  the 
coverage  with  domains  is  drastically  changed  by  annealing. 
So  it  is  likely,  that  during  oxidation  some  defects  at  the 
interface  are  generated,  which  to  some  extent  may  be 
annealed  in  nitrogen. 
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posttreatment  of  the 
Sl(lll)/S102-lnterface 

none 

800  °C 

10  rain 

1000  °c 

10  min 

1000  °C 

60  min 

coverage  of  the 
top-layer  6  C%  of  ML] 

36  ±  6 

34  ±  7 

31  ±  11 

36  ±  5 

asperity  height  A  C&J 

LSI  ±  0.06 

1.48  ±  0.06 

1.44  t  0.16 

1.51  t  0.03 

terrace  correlation 

length  Inm] 

6.5  ±  0.7 

6.5  ±  0.7 

B 

6.5  ±  0.7 

average  terrace  length 
of  the  bottom-layerlnm] 

19  ±  5 

20  ±  6 

22  ±  10 

19  ±  4 

average  terrace  length 
of  the  top-layer  [nml 

10  ±  2 

10  ±  2 

10  ±  3 

10  ±  2 

step  atom  density  1X3 

2.3  ±  0.4 

warn 

2.1  ±  0.6 

2.3  ±  0.4 

~T k  .  A 


posttreatraent  of  the 
Si(lll)/Si02-interface 

none 

800  °C 

10  min 

1000  °C 

10  min 

1000  °c 

60  min 

silicon  coverage  0S1 

IX  of  ML! 

87  ±  2 

76  t  5 

8S  ±  5 

95  ±  2 

domain  correlation 

length  U 

1.16  ±  0.0S 

1.34  ±  0.07 

1.91  ±  0.13 

2.60  ±  0.23 

average  length  of  the 

silicon  domains  [nral 

10  ±  2 

6  ±  2 

14  ±  5 

57  ±  27 

average  length  of  the 

defect  domains  Cnm] 

1.5  ±  0.1 

1.8  ±  0.2 

2.5  ±  0.3 

2.9  t  0.3 

defect  domain 
density  C%] 

2.8  ±  0.5 

4.1  ±  0.8 

u 

■ 

TVfc-Q 


In  earlier  experiments  with  low  resolution  instruments  a 
distinction  between  the  two  types  of  defects  were  not 
possible.  Therefore  the  reported  annealing  effect  on  half 
width  at  the  out-of-phase  condition  (S  =  5.5)  may  be  just 
the  annealing  of  the  defects. 

The  influence  of  the  two  types  of  defects  (steps  and 
inhomogeneities)  on  electrical  properties  so  far  only  can  be 
guessed  /l/.  The  reported  effect  on  mobility  may  be  due  to 
both  types  of  defects.  It  has  been  reported,  that  the 
interface  state  density,  which  smaller  that  the  step  atom 
density  by  orders  of  magnitude,  varies  also  with  the  half 
width.  Here  it  is  more  likely,  that  the  inhomogeneities  are 
the  origin  of  the  sites  for  fixed  charge.  For  a  clear 
conclusion  dedicated  experiments  are  neeaed.  So  far, 
however,  it  is  already  shown,  that  two  different  types  of 
defects  may  be  separated  experimentally,  so  that  a  direct 
correlation  to  electrical  properties  may  be  studied  in 
detail . 


5.  The  Si02/Si(100)  interface 

The  Si (100)  surface  has  for  electron  diffraction  the  disad¬ 
vantage,  that  subsequent  layers  belong  to  different  sublat¬ 
tices  and  are  therefore  not  equivalent.  The  interference 
pattern  may  be  modified  by  different  scattering  factors,  up 
to  the  limit  of  very  different  form  factors  with  a  misinter¬ 
pretation  of  double  steps.  Since  the  two  levels  are 
identical  after  a  rotation  of  90°  ,  the  form  factors  are 
identical,  if  the  incident  beam  is  selectd  with  an 
azimuthal  direction  just  45°  to  the  two  levels.  To  check  the 
influence  of  different  form  factors  onto  the  intensity  for 
the  central  spike,  a  four  level  system  (with  all  levels  at 
same  probabiliy  has  been  calculated  for  a  variing  ratio  of 
the  two  form  factors  (fig.  5.1).  It  is  seen,  that  the 
central  spike  is  not  affected  for  energies  close  to  an  in¬ 
phase-condition.  Therefore  the  analysis  for  rms  roughness 
is  not  affected.  Whereas  at  out-of-phase  for  very  different 
form  factors  (the  ratio  A  of  the  form  factors  A  =  0.0)  the 
high  central  spike  might  look  like  in  a  double  step 
structure,  already  a  factor  of  0.5  brings  the  relative  peak 
intensity  down  to  about  10  %.  It  is  therefore  important  to 
choose  a  diffraction  condition,  so  that  both  levels  contri¬ 
bute  . 

Samples  have  been  oxidized  and  annealed  the  same  way  as  the 
(111)  samples  (10  r.m  oxide  in  dry  oxygen,  increasing 
anneal  in  nitrogen).  A  result  is  shown  in  fig.  5.2  .  It  is 
seen  in  the  upper  half  (direction  for  equivalent  levels) , 
that  four  levels  are  needed  for  description  with  approxi¬ 
mately  the  same  probability  for  each  level.  Close  to  out-of- 
phase  condition  no  central  spike  was  observed,  which  may 


be  also  due  to  a  vanishing  form  factor.  It  is  more 
likely,  however,  that  it  is  due  to  destructive  interference. 
If  the  azimuth  of  the  incident  beam  is  changed  by  45° 
degrees,  the  form  factors  of  the  two  levels  are  quite 
different,  so  that  at  out-of-phase  condition  a  large  central 
spike  is  seen.  Since  that  difference  is  known,  neverthe¬ 
less  the  four  levels  are  to  be  taken  for  description. 

The  evaluation  has  been  done  the  same  way  as  with  the 
Si(lll)  samples  (table  3). 

The  central  spike  provides  the  roughness  (always  1.8  A)  and 
the  number  of  layers  (always  4) .  The  spot  profile  again 
consists  of  a  narrow  Lorentzian  shoulder  (steps)  and  a  broad 
Lorentzian  shoulder  (inhomogeneities)  ,  which  are  given  in 
the  table  with  the  respective  average  widths.  Since  the  data 
have  not  been  collectd  for  a  number  of  energies  as  high  as 
with  the  (111)  samples,  the  interpretation  of  the  broad 
shoulder  as  due  to  inhomogeneities  is  just  done  analogous 
to  the  (111)  surface. 

Since  the  step  height  at  the  (100)  surface  is  only  1.34  8 
(compared  to  3.13  8  at  (111)  ),  the  rms  roughness  of  1.8  8 
over  4  layers  may  yield  a  smoother  surface  on  an  atomic 
scale  than  the  (111)  surface  (with  rms  roughness  of  1.5  8). 
At  the  (100)  interface  the  annealing  both  reduces  the  step 
density  and  the  inhomogeneity  distribution  which  is  in 
contrast  to  the  (111)  results,  reported  above. 

The  present  results  provide  the  first  detailed  description 
of  the  SiOs/SidOO)  interface.  In  TEM  studies  mostly  only 
two  levels  are  seen,  since  the  resolution  is  only  close  to 
1.7  8.  The  present  result  has  the  same  good  interfaces,  it 
resolves,  however,  due  to  the  much  better  resolution,  all 
four  levels,  which  contribute  to  the  formation  of  the 
interface . 


6.  Regular  Step  arrays  at  the  SiOa/Sidll)  interface 

A  set  of  wafers  was  misoriented  away  from  the  dll) 
direction  by  1.7°.  For  a  perfect  interface  a  regular  step 
array  with  a  terrace  width  of  about  100  A  is  expected.  It 
has  been  assumed,  however,  that  the  interface  roughness 
should  obscur  the  regular  array  completely.  Since  some  of 
the  wafers  have  been  carefully  annealed  after  oxidation,  it 
was  of  interest  if  any  regularity  would  be  visible. 

Fig.  6.1  shows  a  clear  splitting  of  the  spot,  which  is 
quantitatively  resolved  in  fig.  6.2  for  a  linear  scan.  The 
clear  splitting  shows  a  good  regularity,  which  is  improved 
by  annealing.  Since  the  same  wafers  have  been  used  for  the 
comparative  study  (see  section  7),  here  only  the  aspect  of 
regular  step  array  is  discussed. 
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Here  for  the  first  time  a  regular  step  array  has  been 
detected  at  the  Si/SiOj  interface.  The  broadening  of  the 
double  spot  and  broad  base  (which  is  important  even  after 
notice  of  the  logarithmic  scale  in  fig.  6.2)  shows 
deviations  from  the  regularity.  For  description  of  those 
irregularities  appropriate  parameters  or  measures  had  to 
be  found.  For  that  purpose  first  a  functional  form  for  the 
terrace  width  distribution  had  to  be  detected.  The  selected 
function  has  to  reproduce  the  correct  spot  profile  not  only 
for  out-of-phase  condition  (fig.  6.2),  the  profiles  at  all 
other  energies  have  to  fit  automatically. 

After  several  attempts  a  special  function,  which  is  derived 
out  of  the  geometric  distribution,  provided  an  excellent 
fit.  The  Gamma  distribution  P(D  =  r*exp(-  r/r)  provides 
the  fit  as  is  shown  in  fig.  6.2.  The  distribution  is  shown 
in  fig.  6.3,  where  the  parameters  a  and  r  have  been 
selected  so  that  in  all  cases  the  average  terrace  width  (and 
therefore  the  splitting  of  fig.  6.2)  is  the  same.  The 
sharper  the  distribution,  the  deeper  is  the  minimum  between 
the  split  spots.  To  reduce  the  number  of  parameters,  only 
two  are  used  for  further  description.  One  is  the  average 
terrace  width,  which  for  a  monotonuous  step  array  is  given 
by  the  angle  of  misorientation.  The  other  one  is  the 
standard  deviation  o  from  the  regular  step  array  with  the 
average  terrace  width.  For  an  ideally  regular  array  it  would 
be  zero.  In  reality  the  terrace  width  varies,  so  that  the 
standard  deviation  is  a  measure  of  irregularity. 


For  the  direction  parallel  of  the  step  edges  the  profile  is 
identical  to  those  found  on  well  oriented  wafers.  Therefore 
the  same  fitting  procedure  and  the  same  parameters  may  be 
used:  the  rms  roughness,  the  average  length  of  an  excursion 
into  the  next  level.  For  the  stepped  surface  the  second 
parameter  may  also  be  expressed  as  the  average  length  of 
straight  portions  (free  of  kinks)  along  the  step  edge 
direction. 

Therefore  the  LEED  evaluation  provides  detailed  information 
on  the  deviations  from  regularity.  All  investigated  samples 
had  the  sane  misorientation,  they  had,  however,  due  to 
different  oxidation  parameters  quite  different  standard 
deviations  from  regularity  and  average  length  of  straight 
edges.  The  results  will  be  discussed  in  section  7  together 
with  the  electrical  and  optical  properties  of  the  same 
wafers . 

Here  only  the  importance  with  respect  to  possible  electrical 
transport  properties  due  to  the  new  periodicity  should  be 
discussed.  Any  new  periodicity  produces  new  Brillouin  zone 
boundaries  and  therefore  new  forbidden  gaps  in  the  band 
structure  (minigaps  and  minibands) .  If  th  mean  free  path  is 
much  larger  than  the  new  period  (10  nm  in  this  case)  the 
carriers  no  more  move  in  the  parabolic  part  of  the 
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conduction  band,  close  to  the  minigaps  drastic  anomalies  in 
conductivity  are  expected,  which  may  be  used  for  a  new 
transistor  concept  (Bloch-transistor) .  The  present  experi¬ 
ments  show,  that  it  is  possible  to  produce  a  regular  step 
array  at  the  Si/Si02  interface.  The  regularity  may  be  varied 
to  a  large  extent  by  appropriate  oxidation  parameters.  It 
would  be  worthwile  to  check  directly  for  those  transport 
anomalies  in  new  experiments. 


7.  Comparative  Study  with  many  methods 
7.1  Samples  and  preparation 

A  set  of  wafers  Si (111)  has  been  oxidized  in  6  groups  with  4 
wafers  each  with  variation  of  parameters,  which  are  shown  in 
the  following  table  7.1. 


Table  7.1 


Group 

I 

II 

III 

IV 

B! 

1 

preoxidation 
(thick  oxide) 

X 

X 

X 

■ 

■ 

- 

oxidation  in  dry  O2 
(thickness  in  nm) 

40 

40 

90 

90 

40 

■ 

wet-dry 

- 

■ 

■ 

■ 

B 

100 

annealing  in  dry  N2 
(1  hr,  1000°C) 

- 

X 

X 

X 

X 

optical 

light  scattering 

before  oxidation 

all  lower  than  6 

after  oxidation 
(%  of  max.  value) 

70 

70 

26 

26 

11 

100 

-  13  - 


The  set  of  parameters  was  selected  to  obtain  informations 
on  the  effect  of  the  following  parameters: 

a)  influence  of  preoxidation  (group  II  and  V) 

b)  influence  of  oxide  thickness  (I  -  IV,  II  -  III) 

c)  influence  of  annealing  in  dry  N2  (I  -  II,  IV  -  III) 

d)  influence  of  wet  oxidation  (VI  -  III) 

The  samples  have  been  selected  with  high  resistivity  and 

polished  on  both  sides,  so  that  ESR  was  possible.  The 

preoxidation  was  on  oxide  of  100  nm  and  an  anneal  in  N2  at 

1000°C  for  1  hr.  This  oxide  has  been  removed  in  HF  before 

the  next  oxidation  was  started.  The  light  scattering  data 
were  recorded  at  Wacker  Chemitronic,  Burghausen.  The  oxide 
thickness  has  been  determined  with  ellipsomtry. 


7.2  Experimental  results 

LEED.  The  wafers  are  those  misoriented  samples,  as  described 
in  section  6,  that  is  a  misorientation  of  1,7°  towards  211. 
As  already  described,  the  regularity  in  the  direction  of  the 
step  array  and  along  the  step  edge  could  be  measured 
independently  by  the  appropriate  scan  direction  of  the 
spot  profile.  The  more  regular  the  step  array  is,  the 
narrower  is  the  probability  distribution  of  the  terraces. 
Its  standard  deviation  from  the  average  terrace  width  is 
therefore  a  measure  for  its  regularity.  For  the  regularity 
along  the  step  edge  the  average  length  of  straight  edges  has 
been  taken.  For  the  samples  of  group  I  the  set  of  data  was 
not  complete  for  a  full  evaluation,  the  data,  however,  were 
sufficient  to  give  the  rank  of  the  quality. 

The  results  are  shown  in  table  7.2. 


Table  7.2 


Group 

I 

II 

hi 

IV 

n 

VI 

standard  deviation  (nm) 

■ 

from  regular  step  array 

- 

1,6 

1,0 

1,4 

1,8 

1 

rank  fo  quality 

6 

4 

1 

3 

5 

H 

average  length  of 
straight  edged 

12 

18 

21 

1 

20 

33 

rank  of  quality 

6 

4 

2 

3 

1 

Electron  spin  resonance 
(spin  /  cm2 ) 

2X101  1 

5x10* o 

- 

- 

- 

Electrical  breakdown 
(MV  /  cm) 

8,6 

9,3 

8,0 

8,6 

9,1 

8,7 

rank 

4 

i 

6 

4 

2 

3 

The  ESR  measurements  have  been  kindly  made  by  E.  Poindexter 
ET  &  DL,  Fort  Monmouth  N.  J.  so  far  for  the  samples  of  the 
group  I  and  XI.  The  other  measurements  are  in  progress.  His 
comments  are,  that  the  spin  density  is  much  lower  than 
expected.  An  explanation  cannot  be  given  so  far.  The 
decrease  from  group  I  to  II  (annealing)  is  anyway  very 
drastic.  Electrical  breakdown  has  been  measured  at  Wacker- 
Chemitronic  for  two  wafers  out  of  each  group.  A  typical 
example  is  shown  for  the  wafers  out  of  group  II  in  fig.  7.1. 
The  average  breakdown  field  strength  and  a  rank  with  respect 
to  this  quality  is  also  shown  in  table  7.2. 


7.3  Discussion 

The  summarized  data  of  table  7.2  show  a  lot  of  interesting 
features:  The  wet-dry  oxidation  is  superior  in  most  proper¬ 
ties,  only  with  the  breakdown  it  is  not  the  best,  here, 
however  it  had  the  lowest  number  of  early  failures.  The 
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other  surprising  feature  is,  that  the  stray  light  inten¬ 
sity  after  oxidation  does  not  correlate  with  the  other 

properties,  which  cannot  be  explained  so  far.  The  other 

parameters,  however,  may  be  discussed  together. 

The  influence  of  preoxidation  (group  II  -  V)  is  fairly 

small,  it  is  not  clear,  if  a  systematic  difference  in 

quality  is  derived  at  all.  Only  the  stray  light  has  a 
strong  effect. 

The  increased  oxide  thickness  (group  I  -  IV,  II  -  III) 
increases  the  regularity  and  decreases  the  stray  light.  Only 
the  extremely  small  breakdown  of  group  III  does  not  fit  into 
this  pattern. 

The  influence  of  annealing  in  dry  N2  at  1000°C  after 
oxidation  gives  a  very  clear  result  (group  I  -  II,  IV 
III) .  Both  the  structural  homogeneity,  the  ESR  spin  density 
and  the  electrical  breakdown  show  an  improvement.  Again  the 
breakdown  of  group  III  is  exceptional.  Here  the  stray  light 
is  not  affected  at  all. 

As  a  surprise  the  wet-and-dry  oxidation  had  a  superior 
performance  with  respect  to  homogeneity.  Although  the 
stray  light  was  high,  th  homogeneity  is  better  than  group 
II,  III  and  V,  which  are  closest  in  preparation  due  to 
annealing . 

The  present  results  for  the  first  time  combine  quite  a 
variety  of  structural  and  electrical  methods  for  a  study  of 
defects  at  the  Si/SiOj  interface.  Due  to  the  novel  evalua¬ 
tion  procedures  of  high  resolution  LEED  measurements  the 
smoothness  is  characterizd  in  much  more  detail  as  before. 
The  stray  light  measurements  were  not  easily  correlated 
here.  It  should  be  pointed  out,  however,  that  recent 
measurements  of  stray  light  with  well  defined  angular 
variations  are  exactly  correlated  with  LEED  data,  so  that 
the  roughness  at  the  interface  is  described  by  the  same 
geometric  distribution  from  atonic  distances  up  to  several 
microns  (publication  in  preparation,  thesis  completed) .  The 
electrical  measurements  include  so  far  only  ESR  and  break¬ 
down  field  strength,  which  already  show  very  relevant 
correlations.  Since  for  the  present  series  high  resistivity 
silicon  has  been  selected  in  favor  of  ESR,  it  was  for  the 
same  reason  not  possible  to  do  CV  measurements  or  radiation 
damage  studies.  Therefore  a  new  series  with  lower  resisti¬ 
vity  has  been  started. 


16 


8.  Preparation  of  a  study  of  a  possible  correlation  of 
defects  and  radiation  damage 

A  set  of  wafers  have  been  prepared  and  oxidized  with 
different  oxidation  parameters,  so  that  a  broad  range  of 
steps  and  other  defects  are  expected.  Since  the  oxidation 
was  not  done  in  our  lab,  we  had  to  accept  long  delays  to  get 
a  high  quality  oxidation  at  the  level  of  industrial 
standard.  So  we  ran  out  of  time  to  do  first  the  LEED 
measurements,  which  would  be  the  prerequisite  for  radiation 
damage  experiments  by  Dr.  B.  Me  Lean  at  the  Harry  Diamond 
Lab.  As  soon  as  we  have  funding,  we  will  continue  this 
series  of  measurements. 


9.  Summary  and  outlook 

During  the  investigation  the  evaluation  procedures  for  Spot 
Profile  Analysis  of  LEED  (SPA-LEED)  has  been  developed  to  an 
unexpected  high  level,  since  due  to  the  now  available  high 
resolution  LEED  systems  qualitatively  and  quantitatively 
a  big  progress  was  possible.  A  surface  is  characterized  now 
laterally  up  to  more  than  100  nm,  completely  new  is  the 
vertical  analysis  providing  the  rms  roughness  and  the 
determination  of  the  levels  contributing  to  the  interface. 
By  a  careful  study  of  energy  dependance  also  a  second  type 
of  defect  (other  than  steps)  could  be  identified  and 
separated . 

With  this  evaluation  procedure  a  precision  analysis  has 
been  made  for  the  (111)  interface,  showing  in  detail  the 
effect  of  annealing. 

Also  the  (100)  interface  now  could  be  analyzed,  to  that  also 
the  technologically  more  relevant  interface  is  now  studied, 
showing  in  principle  the  same  dependance. 

The  first  big  comparative  study  included  besides  LEED  also 
ESR,  optical  stray  light  and  electrical  breakdown,  showing 
the  importance  of  the  defect  concentration  (as  determined  by 
LEED)  for  electrical  properties.  The  result  suggest  to  look 
for  more  correlations  between  defect  concentrations  and 
other  important  electrical  properties  like  radiation  hard¬ 
ness  or  mobility  in  the  inversion  layer,  which  should  be 
done  in  a  subsequent  study. 
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Figure  Captions 

Fig.  4.1a)Spot  Profile  for  a  SiOj/Si(lll)  interface  for  out-of¬ 
phase  condition.  The  solid  line  corresponds  to  the 
expected  profile  of  a  perfect  crystal  (instrumental 
broadening) 

b)same  as  a)  with  the  two  Lorentzian  shoulders  shown  as 
solid  lines 

Fig.  4.2a)Same  as  fig.  4.1  for  in-phase  condition.  The  existence 
of  a  broad  shoulder  indicates  the  presence  of  defect 
areas  (other  than  steps) 

b) same  as  a)  with  the  broad  Lorentzian  shoulder  as 
a  solid  line 


Fig.  4 . 3a) Fraction  of  the  integrated  narrow  shoulder  Is  out  of 
the  total  intensity  VS  energy, 
b)same  as  a)  for  the  broad  shoulder 


Fig.  4.4  Fraction  of  the  integrated  central  spike  intensity 
Ip  out  of  the  sum  of  Ip  and  Is  (=  narrow  shoulder) , 
to  determine  the  vertical  structure 


Fig.  5.1  Fraction  of  the  integrated  central  spike  out  of  the 
total  intensity,  calculated  for  a  Si (100)  as  a 
4  level  system  with  same  probability  in  each  level. 
A  is  ratio  of  the  form  factors  of  the  two  non 
equivalent  levels. 


Fig.  5.2  Measured  fraction  of  central  spike  for  a  Si (100) 

interface  with  an  azimuthal  incident  direction,  so 
that  A  =  1  (upper  half)  or  A  s  0  (lower  half) 


Fig.  6.1  Area  scan  of  a  split  diffraction  spot  due  to  a  nearly 
regular  step  array  for  an  energy  close  to  an  out-of¬ 
phase  condition  (both  peaks  nearly  equal) 


Fig.  6.2  Linear  scan  through  the  two  peaks  (with  integration 

over  the  direction  perpendicular  to  it)  of  fig.  6.1 
for  out-of-phase  condition. 


Fig.  6.3  Calculated  terrace  probability  distribution  for  the 
Gamma-function  with  always  the  same  average  terrace 
width  (T  =  28)  and  different  standard  deviation. 


Fig.  7.1 


Breakdown  of  oxides  on  wafers  of  group  II 
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